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CHAPTER  I 


DEVELOPMENT  OF  PULSED  RADAR  SIGNAL  PROCESSING 

The  Theory  of  Matched  Filter  Signal  Processing 
The  basi„  concept  of  matched  filtering  evolved  from  an  effort  to 
obtain  a  better  theoretical  understanding  of  the  factors  leading  to 
optimum  performance  of  radar  systems.  Matched  filtering  is  the  optimum 
linear  processing  for  radar  signals  (Cook  and  Bernfeld  1967).  It 
transforms  the  raw  radar  data,  at  the  receiver  input  and  in  the  presence 
of  White  Gaussian  noise,  into  a  form  that  yields  optimum  detection 
decisions.  Some  of  these  decisions  might  be  target  or  no  target  or 
for  estimating  target  parameters  of  range  and  velocity  with  minimum  rms 
errors,  and  also  for  obtaining  maximum  resolution  among  a  group  of 
targets . 

The  characteristics  of  matched  filters  can  be  designated  by  either 
a  frequency  response  function  or  a  time  response  function,  each  being 
related  to  the  other  by  a  Fourier  transform  operation.  In  the  frequency 
domain,  the  matched  filter  has  a  transfer  function,  H(u)) ,  that  is  the 
complex  conjugate  of  the  spectrum  of  the  signal  that  is  to  be  processed. 

-JwTd 

H(u.')  =  KS*(w)e  (1) 

where  S(w)  is  the  spectrum  of  the  input  (received)  signal,  s(t),  and 
T^  is  a  delay  constant  required  to  make  the  filter  physically  realizable. 
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The  K  is  a  normalizing  factor  required  to  yield  a  unity  gain.  K  and 
can  be  ignored  in  formulating  the  underlying  relationships. 

H(oi)  =  S*U) 

The  corresponding  time  domain  relationship  is  obtained  from  the 
inverse  Fourier  transform  of  H(u).  The  result  being  that  the  impulse 
response  of  the  matched  filter  is  the  time  inverse  of  the  known  signal 
function.  If  h(t)  represents  the  matched  filter  impulse  response, 
the  relationship  can  be  written  as  follows: 

h(t)  =  KS(Td  -  t)  (2) 

Again,  ignoring  K  and  Trf,  the  basic  relationship  is: 

h(t)  =  S(-t)  (3) 

This  matched  filter  relationship  can  be  derived  through  optimization 
of  several  different  characteristics.  One  method  of  deriving  the 
matched  filter  relationship  is  to  optimize  the  signal-to-noise  ratio. 

The  following  is  a  reiteration  of  this  development  as  shown  by  Cook 
and  Bernfel d  (1967). 

Tee  following  section  is  devoted  to  showing  that  the  matched  filter 
is  the  optimum  signal  processing  to  maximize  the  signal-to-noise  ratio. 
The  signal-to-noise  ratio  ( S/N )  can  be  defined  as: 

1  r-  Le.iL  Instantaneous  Output  Signal  Power 

M  Hi  if  mit  Ma  i  c  o  Dnuftr* 

It  uUl'k’UV  1  I 


(4) 
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The  objective  here  is  to  determine  the  properties  of  the  linear  system 
(filter)  that  maximize  the  S/N  ratio  as  defined  previously.  By 
definition,  the  peak  instantaneous  output  signal  power  is  the  square 
of  the  maximum  signal  output  voltage,  g(t),  in  the  absence  of  noise 
and  assuming  an  output  impedance  of  1  ohm. 


g(t)  =  /  S(u)  H(o»)  ejVt  dco  (5) 

-oo 

where  S(<jj)  is  the  signal  spectrum  and  H(w)  is  the  system  transfer 
function.  Let  t  =  be  the  time  when  g(t)  is  desired  to  be  a  maximum, 
where  is  equal  to  or  greater  than  the  signal  duration.  Then: 

g(Td)  =  jb  /  S(o>)  H(w)  e  daj  (6) 

— oo 

The  normalized  noise  power  at  the  filter  output  is  expressed  as: 


a 


2 


N  ®  , 

5s-  /  |H(«)|2 

—  OO 


(7) 


where  NQ  is  the  single-sided  noise  power  density,  in  watts/cycle/second, 
at  the  filter  input.  The  S/N  ratio  to  be  maximized  is: 


S 

N 


g2  <V 


2 

a 


(8) 


The  optimization  will  be  accomplished  in  terms  of  finding  the  best 


H(ui]  when  all  other  parameters  are  fixed.  This  can  be  accomplished 
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by  using  Schwarz's  inequality.  The  signal-to-noise  ratio  to  be 
maximized  is  then: 


S 

N 


.  «  jwTj  i 2 

2^  f  S(u)  H(w)  e  dw 

—CO  _  _ 

N  ®  ? 

/  1H{m)  |2  do, 


(9) 


where  Schwarz's  inequality  is  given  by: 


(10) 


S  >  xl<“> 

-00 


X2(w)  du 


<  7  |Xj(u)  I2  du  ±  /|X2(«i)|2  du) 


where  X^w)  and  X2(w)  are  two  arbitrary  functions. 
Equation  (10)  can  be  rewritten  as: 


/  X^ (u) )  X2(w)  dw 


/| Xj(w)  |  du 


OD  2 

/  |X2(u)l  dw 

-CO 


(ID 


The  right  side  of  equation  (11)  is  maximized  when  (Gradshteyn  and 
Ryzhik  1980): 


If: 


Xj(u)  =  kj  X*(u) 


(12) 


Xj (u) )  =  ^  |H(u)|2 


(13) 


and: 


X2 (oj)  =  S (w)  H(ai)  e 
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Therefore,  applying  the  complex  conjugate  and  time  delay 
properties  of  the  Fourier  transform,  the  impulse  response  of 
the  matched  filter  is  given  by: 

h{t)  =<f_1  (HU)}  -  K  $(Td  -  t)  (19) 

Pulse-Compression  Concepts 

The  concept  of  pulse-compression  was  conceived  out  of  a  need  to 
improve  on  the  performance  of  pulsed  radar  systems.  Near  the  end  of 
World  War  II,  the  need  for  radically  improved  radar  performance  became 
apparent  (Cook  1960).  It  also  became  obvious  that  Gated  RF  could  not 
provide  this  improved  performance  due  to  power  limitations  of  the 
transmitters.  These  limitations  existed  in  two  areas.  The  peak  power 
available  in  the  transmitter  tubes  was  limited,  and  secondly,  many  of 
the  transmitter  components  could  not  have  operated  under  any  higher 
power  conditions.  One  possible  solution  to  the  problem  would  have  been 
to  increase  the  average  power  capabilities  of  the  transmitter  tubes 
by  going  to  wider  pulse  widths.  This  would  provide  greater  range 
capabilities,  but  the  conflicting  requirement  of  improved  resolution 
for  ground  mapping  and  separation  of  aircraft  targets  in  massed 
formations  would  be  undermined.  Therefore,  this  most  obvious  method 
of  improvement  was  not  acceptable. 

The  proposed  solution  to  the  problem  of  trying  to  improve  radar 
performance  was  conceived  by  several  individuals  serving  on  both  sides 
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during  the  war  (Cook  and  Bernfeld  1967).  The  solution  basically 
provided  that  a  wide  pulse  be  transmitted,  during  which  the  carrier 
frequency  would  be  linearly  swept.  This  was  seen  to  yield  a  correlation 
between  time  and  frequency  that  could  be  exploited  in  the  receiver. 

Tiiis  would  be  accomplished  by  using  a  filter  having  a  linear  time 
delay  versus  frequency  characteristics  such  that  it  would  apparently 
delay  one  end  of  the  received  wide  pulse  by  a  greater  amount  than  the 
other,  thus  causing  the  signal  to  compress  in  time  and  increase  in 
peak  amplitude.  Thus  was  born  the  concept  of  pulse-compression.  Actual 
implementation  did  not  occur  until  the  early  1950s  when  the  high  power 
Klystron  was  developed. 

Gated  RF  Pulsed  Radar  Techniques 
The  concept  of  the  gated  RF  pulse  was  dominant  in  radar 
applications  until  the  post-World  War  II  period  (Barton  1976).  The 
basic  concept  is  to  transmit  a  short  burst  of  high  energy,  uniform 
amplitude  and  constant  frequency  RF  (radio  frequency)  signal  and  wait 
for  a  return  signal  reflected  off  a  target.  A  matched  filter  is  used 
in  the  receiver  to  process  the  return  signal  and  yield  an  output 
signal  that  contains  valuable  information  concerning  target  range  and 
velocity.  A  block  diagram  of  the  receiver  matched  filter  is  shown  in 
Figure  1. 
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G(w)  =  the  spectrum  of  the  output  signal 
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The  output  auto-correlation  function,  g(t),  for  gated  RF  can 
be  obtained  from  a  very  simple  and  straightforward  closed  form 
analysis.  The  received  signal,  s(t),  shown  in  Figure  2,  can  be 
expressed: 


s(t)  =  cos  ui  t  0  <  t  <  T 
o  —  — 


This  can  be  written  in  terms  of  a  rect  function  as: 


t  - 


s(t)  =  cos  w  t  •  rect  ( — 
o  T 


Now,  applying  simple  Fourier  transform  pairs  yields  an  expression  for 
the  spectrum,  S(f),  as  follows: 


S(f)  =  Vi  6(f  -  f0)  +  >s  5 ( f  +  fQ)]  *  [T  SINC  fT] 

where  *  denotes  the  convolution  in  the  frequency  domain.  Completing 
the  convolution  will  yield  frequency  shifted  SINC  functions  as  follows: 

S(f)  =  \  SINC  (f  -  f0)T  +  \  SINC  (f  +  f0)T  (20) 

Here,  it  will  be  assumed  that  fQ  is  a  great  enough  value  such  that  the 
negative  frequency  portion  has  negligible  effect  in  the  positive  frequency 
range.  Figure  3  shows  a  plot  of  the  positive  spectrum.  The  spectrum 

r  t  MV 

is  the  familiar  — jj—  shape  with  -13  dB  side  lobes.  The  center  frequency 


vyyyy 
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(carrier),  f  ,  is  70  MHz.  This  same  carrier  frequency  will  be  used 
throughout  this  paper. 

Since  all  systems  considered  in  this  paper  are  matched  filter 
systems,  the  matched  filter  case  will  be  investigated  here.  The  matched 
filter  impulse  response,  h(t),  can  be  written  as: 


Now  with  simple  trigonometric  identities,  h(t)  becomes: 

t-J 

h(t)  =  [cos ( 2rrK)  cos((jjQt)  +  sin(2rK)  sin (co0t ) ]  •  rect  ( — j-=-) 
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Since  cos(2ttK)  =  1  and  sin(27iK)  =  0,  this  becomes: 

h { t )  =  cos(to0t)  *  rect  ( — ^-=-) 

This  is  identical  to  the  received  signal,  s(t).  Therefore,  the  matched 
filter  spectrum,  H(u)),  will  be  identical  to  S(w).  The  reader  is 
referred  to  figures  2  and  3  for  plots  of  these  functions. 

The  desired  output  from  the  matched  filter,  the  auto-correlation 
function,  can  be  obtained  analytically  by  doing  the  actual  convolution 
where: 


g(t)  =  s(t)  *  h(t) 

oo 

-I  s(t)  h(t  -  t)  dT 

— oo 

or  by  finding  G(ui)  and  taking  the  inverse  Fourier  transform: 

6(f)  =  [S(f)]? 

which,  considering  only  positive  frequencies,  can  be  written  as: 

t2  o 

G(f )  =  —  SI  No  [(f  -  fc)T] 

Now  using  a  very  common  Fourier  transform  pair,  g(t)  can  be  written 
directly  as: 

r  £  cos  ui  t  0  <  t  <  T 

r,  o 

g(t)  =  j  (T  -  |-)  cos  wQt  T  <_  t  <  2T 


0 


otherwise 


13 


Figure  4  is  a  plot  of  the  magnitude  of  g ( t)  obtained  from  the  Surface 
Acoustic  Wave  Computer-Aided  Design  program  (SAWCAD)  by  squaring  the 
signal  spectrum,  S(f)  (see  Figure  3)  and  then  taking  the  inverse  FFT 
(Fast  Fourier  Transform).  Due  to  the  FFT  processing,  only  the 
positive  portion  of  the  signal  is  plotted. 
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Figure  4.  Auto-Correlation  Function  for  a  Gated  RF  Pulse  Matched 
Filter  System  Output. 


Linear  FM  Pulse-Compression  Systems 
Initial  implementation  of  pulse-compression  was  in  the  form  of 
linear  FM  (Klauder  et  al.  1960).  It  is  currently  the  most  widely  used 
type  of  pulse-compression  waveforms.  Linear  FM  pulse-compression  is 
easy  to  implement  and  for  a  majority  of  radar  applications  it  has  the 
best  overall  properties.  In  general,  linear  FM  waveforms  provide 
improvement  over  Gated  RF  by  satisfying  five  basic  requirements. 

These  requirements  are  (Ramp  and  Wingrove  1961): 

1.  A  product  of  bandwidth  and  duration  greater  than  unity 
(for  good  resolution  of  both  range  and  velocity) 

2.  A  rectangular  amplitude  spectrum  (for  best  range 
resolution) 

3.  A  rectangular  envelope  (for  best  velocity  resolution) 

4.  Phase  dispersion  (that  is,  a  non-linear  phase  spectrum 
to  permit  pulse-compression) 

5.  Phase  matching  can  be  accomplished  using  a  single  filter 
Here,  linear  FM  pulse-compression  will  be  examined  in  the  same 

type  of  matched  filter  system  as  the  Gated  RF  was.  The  transmitted 
waveform  for  the  linear  FM  pulse-compression  signal  is  written: 

s(t)  =  cos to>0 1  +  ^  t2]  “  \  <  t  <  ■£  (21 

=  0  elsewhere 

where : 

u0  =  carrier  frequency  (radians/sec) 
v  =  rate  of  frequency  sweep  (Aw/T) 

Aw  =  swept  spectrum  bandwidth  (yT) 

T  =  uncompressed  pulse  width 
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Figure  5  is  a  plot  of  a  sample  linear  FM  signal. 

The  spectrum  of  the  linear  FM  signal  can  be  expressed  by  the 
fol lowing: 

T 


S(u>) 


COS[(i)0t  + 


li^-l  e'j”1 


dt 


T  .2 

2  j[(o)n  +  (22) 

=  k  /  e  0  1  dt 

T 

'  7 


1  t2 

2  -j[U.  *u)t  + 

+  /  e  0  £  dt 

T 
2 


where  the  first  integral  defines  the  positive  frequency  spectrum  and 
the  second  integral  defines  the  negative  frequency  spectrum.  Since  the 
negative  frequency  spectrum  has  a  negligible  contribution  at  positive 
frequencies,  it  will  be  ignored  here.  By  completing  the  square  of 
the  bracketed  term  in  equation  (22),  the  spectrum  becomes: 


l-M- 


(u)  -  oj  )‘ 


S(u>.)  =  H  e 


2u 


-)1  2 
/ 


y 

■ST 

C 


(+  - 


O)  -  to 


°)2 


dt 


(23) 


Using  the  change  of  variable: 
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Figure  5.  Plot  of  a  Sample  Linear  FM  Transmitted  Signal. 


03  -  to 

fii  (t - - — °)  =  Sn  X 


and 


dt  -  -  dx 

V 


The  spectrum  can  now  be  written: 


S  (to )  =  *5  rr  e 


p 

(u)  ”  U)  )  y  Y 

ZE-®-*  2  j 


~2~ 

f  e  c  dx 
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where : 


y  _  +  (U)  ~  t°o) 

1 


and 


X  =  ^  ^  '  ajQ) 

2  v'v'u 


The  spectrum  can  now  be  expressed  in  terms  of  the  Fresnel  integrals 
as : 

(w  -  u  )2 

-Jt-nrr2-) 


S(w)  *  H  ’  ~  e 


2u  ;  [C{Xj)  +  jS(Xx)  +  C(X2)  +  jS(X2)l 


(25) 


The  above  expression  for  the  linear  FM  spectrum  can  be  broken  down 
into  three  components: 

Amplitude  Term: 


|S(u>)j  =  h  ^  UC(X1)  +  C(X2)]2  +  IS(XX)  +  s(x2)]2)Js 


Square  Law  Phase  Term: 


(»0  -  «)‘ 


1  2p 


and : 
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Residua]  Phase  Term: 

_J  S(XX)  +  S(x2) 

^2  =  “  tan  [C(XX)  +  C(k2')] 

where : 

X  2 

C(X)  =  /  cos  dy 

o  c 

and 

X  2 

S(X)  =  /  sin  dy 

o  c 

are  the  well-known  Fresnel  integrals. 

An  illustrated  plot  of  the  received  signal  spectrum,  S(oj),  is 
shown  in  Figure  6.  The  actual  plotting  was  accomplished  using  the 
FFT  and  graphics  capabilities  of  SAWCAD. 

Now,  using  matched  filter  principles,  the  impulse  response  of 
the  matched  filter  can  be  written  directly  as  the  time  inverse  of  the 
transmitted  signal. 

h(t)  =  ^  ^  cos[u>ot  -  -  \  <  t  £  Y  (26) 

=  0  elsewhere 

where  /2y/iT  is  the  normalizing  factor  required  to  give  the  filter  unity 
gain.  An  illustration  of  the  matched  filter  impulse  response,  h(t),  is 
shown  in  Figure  7. 


& 
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To  obtain  the  matched  filter  frequency  spectrum,  the  matched 
filter  relationship,  H(w)  =  S*(u>)>  is  used  to  write  H(u)  as: 


_j($i  +  <j>2^ 

H(to)  =  S*(co)  =  |  S  (to)  |  e 


/  TT 


J  {[C(X1)  +  C(X2)]l  +  [S(Xj)  +  S(X2)n  e 


-j(- 


{(jj-w  )‘ 


2p 


(27) 


where  <j>j  is  the  square  law  phase  term: 

(“o  •  “>* 
*1  "  2y 


(28) 


and  <j>2  is  the  residual  phase  term: 

S(Xj)  +  S(X2) 
^2  5  "  tan  ^(Xj)  +  C(X2)] 


(29) 


The  residual  phase  term  can  be  neglected  here  because  it  can  be 
approximated  by  a  constant  phase  term  over  the  spectrum  bandwidth,  Af. 
Figure  8  is  an  illustration  of  the  matched  filter  spectrum,  H(co). 

The  output  of  the  matched  filter  is  the  auto-correlation  function 
of  the  input  signal,  s(t).  This  function  can  be  obtained  by  convolving 
the  input  signal,  s(t),  with  the  matched  filter  impulse  response, 
h(t) ,  yielding: 


<jo 

g(t)  =  /  S(t)  h(T  -  t)  dt 


\  '  *, 


P 


£ 

A 

A' 


M- 

A 


m 


L 


i 


§ 


>** 


-00 


(30) 
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Figure  8.  Linear  FM  Matched  Filter  Frequency  Response. 
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Another  method  is  to  multiply  the  input  signal  spectrum  by  the  matched 
filter  spectrum  and  then  take  the  inverse  Fourier  transform,  yielding: 


g(t)  =  V1  {G(tj)J  =  V1  {s(u)  ■  H(u>)>  (31) 

Using  this  second  method,  g(t)  can  be  derived  in  the  following  manner: 

g(t)  =  0-1  (S(u>)  *  H(w)}  =  r1  {6(w)>  (32) 


(32) 
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where: 


G(u>)  =  S(uj)  ’  H(w)  =  S(w)  •  S*(gj)  (33) 

■  |S(U)|2 

Substituting  for  S(io)  yields: 

G(u)  =  \h  {[C(Xj)  +  C(X2)]2  +  [S(Xj)  +  S(X2)]2}Js|2  (34) 


fcJUCtXj)  +  C(X2)]2  +  [S(X1)  +  S(X2)]2} 


(35) 


Taking  the  inverse  Fourier  transform  yields  the  following  result: 

sin  (T  ~  1*1)3 


g(t) 


uT 


cos  (to  t)  -T<  t<  T  (36) 


This  equation  now  represents  the  usable  output  from  the  matched 
filter  and  with  further  receiver  processing  can  provide  information 
concerning  range,  velocity  and  dimension  about  a  target.  By  definition, 
Af  =  pT/2n  and  1/TAf,  where  the  latter  is  the  amount  the  output  pulse 
width  has  been  reduced.  The  factor,  TAf,  defines  a  pulse  width  ratio 
of  the  input  and  output  functions  known  as  the  compression-ratio.  This 
same  factor  represents  the  amount  the  peak  power  level  has  increased 
over  the  Gated  RF  case.  Figure  9  is  a  plot  of  the  auto-correlation 
function  for  a  sample  linear  FM  signal. 
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TIfE  in  microSECOWft 


Figure  9.  Auto-Correlation  Function  for  a  Matched  Filter  Linear  FM 
Pulse-Compression  System. 


Linear  FM  pulse  compression  has  provided  the  necessary 
improvement,  in  range,  velocity  and  resolution,  over  Gated  RF.  But, 
linear  FM  still  has  its  shortcomings.  One  such  weakness  is  that  the 
auto-correlation  function  has  fairly  high  sidelobes;  these  sidelobes 
increase  as  the  compression  ratio  is  increased.  Several  filtering 
techniques  have  been  employed  to  reduce  these  sidelobes.  One  method 
used  to  suppress  these  sidelobes  is  to  deliberately  mismatch  the 
amplitude  response  of  the  matched  filter  (Bernfeld  et  al .  1965).  This 
has  the  effect  of  lowering  sidelobe  levels,  but  introduces  a  1-2  dB 
mismatch  loss. 
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Non-Linear  FM  Pulse-Compression  Techniques 
In  an  effort  to  overcome  the  mismatch  loss,  but  still  achieve  low 
auto-correlation  si  delobe  levels,  non-linear  FM  has  been  investigated. 
Non-linear  FM  functions  with  odd  symmetry  have  been  investigated 
(Bernfeld  et.  al .  1965),  which  yield  lew  sidelobe  levels  while  avoiding 
the  mismatch  loss.  One  such  class  of  functions  yields  spectra  at  the 
mismatch  filter  output  given  by: 

|G(w)  |=  cos  [ — ~ - ]  for  |w0-w|<^ 


The  matched  filter  auto-correlation  waveform  of  these  signals  is: 


9(t)  = 


1_ 
2t r 


,,  j.  Au) 

“o  +  ‘7 

/  cosn  [ - 2- 

«0  -  ^  Au) 


-  ui)  j(u>rt  -  to)  td  (u)0  -  0)) 


-1  e 


The  authors'  results  are  given  in  Table  1  (Cook  et  al .  19  ). 


TABLE  1 

AUTO-CORRELATION  SIDELOBE  COMPARISONS  FOR  NON-LINEAR 
FM  WITH  COSINE"  OUTPUT  SPECTRA 


n 

Aw 

n 

Aojj. 

SIDELOBE 

1 

1.00 

-23  dB 

2 

1.21 

-32  dB 

3 

1.40 

-40  dB 
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It  should  be  noted  that  non-linear  FM  is  only  suitable  for  applications 
in  which  the  expected  range  of  Doppler  frequency  shifts  is  quite  small. 

Another  example  of  using  a  non-linear  FM  waveform  to  achieve  low 
time  sidelobes  has  been  presented  in  which  the  author  was  able  to 
achieve  -27  dB  sidelobes  in  a  system  with  a  time-bandwidth  product  of 
22  and  a  0.1  dB  mismatch  loss  (Mill ett  1970).  This  was  achieved  by 
choosing  a  particular  weighting  function  of  the  expanded  pulse  spectrum. 
In  particular: 


|F(«)|»  4  4  (l  -  K)  cos2 

where  K  is  the  amplitude  of  the  output  spectrum  and  A  is  the  signal 
bandwidth.  Then,  by  employing  a  receiver  with  an  amplitude  response 
that  has  the  same  variation  with  frequency  and  a  dispersion  of  the 
opposite  sense  as  that  in  the  received  signal,  the  system  is  a  matched 
filter  system  with  "cosine  squared  on  a  pedestal"  weighting.  Each 
dispersive  filter  in  this  case  was  synthesized  as  a  cascade  of 
bridged-tee  networks. 

Both  of  these  techniques  proved  to  be  usable;  however,  the 
design  methods  are  complicated  and  iterative.  In  the  next  section, 
a  new  non-linear  FM  pulse-compression  waveform  is  investigated  for 
design  on  a  surface  acoustic  wave  device. 


CHAPTER  II 


OBJECTIVE  OF  PROPOSED  WORK 

The  main  objective  of  this  thesis  is  to  introduce,  analyze  and 
evaluate  a  new  non-linear  FM  pulse-compression  technique  to  be 
implemented  on  a  surface  acoustic  wave  (SAW)  device.  Although  the 
implementation  has  not  actually  been  accomplished,  it  is  intended 
that  it  be  done  using  a  computer-aided  design  program,  SAWCAD, 
developed  at  the  University  of  Central  Florida.  Design  procedures 
and  considerations  will  be  discussed  (Malocha  and  Richie  1984). 
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CHAPTER  III 


INTRODUCTION  AND  EVALUATION  OF  NEW  NON-LINEAR 
PULSE-COMPRESSION  TECHNIQUE 

Theory  and  Closed  Form  Analysis 

The  principle  benefit  of  using  a  non-linear  modulating  term  is 
the  reduction  in  auto-correlation  sidelobe  level  and  a  like  gain  in 
mainlobe  signal  strength  in  a  matched  filter  pulse-compression  system. 
It  has  been  shown  that  for  many  applications  in  which  the  range  of 
Doppler  shifts  is  quite  small  (stationary  target),  the  signal  auto¬ 
correlation  waveform  is  a  major  consideration  (Bernfeld  et  al.  1965). 
The  desired  auto-correlation  waveform  should  exhibit  very  low  si  delobes 
in  order  to  reduce  the  masking  effect  of  large  signals  on  small  signals 
net  at  the  same  range.  To  obtain  low  sidelobes  with  linear  FM  would 
require  using  a  frequency  weighting  mismatch  filter  after  the  matched 
filter.  This,  in  turn,  leads  to  a  reduction  in  signal-to-noise  ratio 
relative  to  the  ideal  matched  filter  case.  The  use  of  a  non-linear  FM 
function  as  a  means  of  obtaining  matched  filter  waveforms  with  low 
range  sidelobes,  while  avoiding  mismatch  loss,  has  been  investigated 
(Cook  et  al .  1965) . 

In  this  section,  a  non-linear  FM  function  is  investigated  for  its 
potential  as  a  low  sidelobe  level  matched  filter  pulse-compression 
system.  This  system  will  again  be  considered  in  the  general  receiver 
configuration,  as  shown  in  Figure  10. 
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Figure  10.  Block  Diagram  of  the  Receiver  Matched  Filter  for  a  New 
Non-Linear  FM  Pulse-Compression  System. 

where: 

s ( t )  =  the  ideal  received  signal  (the  transmitted  signal) 

S(oj)  =  the  spectrum  of  the  received  signal 
h(t)  =  the  matched  filter  impulse  response 
H(oj)  =  the  matched  filter  spectrum 
g(t)  =  the  auto-correlation  function 
G(oj)  =  the  output  spectrum 

Here,  the  non-linear  FM  waveform,  s(t),  is  obtained  by  using  a 
form  of  the  class  of  the  raised-cosine  family  of  functions  known  as 
the  Eigen  function.  In  particular,  the  hj^t)  Eigen  function  has 
been  utilized,  where: 
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This  function  is  shown  in  Figure  11,  plotted  from  -  ^  to  Due  to 
the  symmetry  about  zero,  only  the  portion  from  -  to  0  is  used  to 
create  the  up-chirped  (expanded  pulse)  received  signal,  s(t).  This 
portion  of  the  Eigen  function  is  shown  in  Figure  12. 

In  order  to  obtain  a  modulated  frequency  deviation,  Af,  that  is 
proportional  to  the  Eigen  function,  it  is  necessary  to  use  the  integral 
of  h3j(t)  as  the  phase  term  in  s{t).  The  non  linear  FM  waveform  is 
then  written  as: 

s(t)  -  cos{(i)Qt  +  Df  [.43t  +  sin(^)  +  sin(^)]} 

for  -  ^  1  t  <_  j  where  T  -  2x 

Now,  utilizing  only  that  portion  of  the  signal  from  -  j  to  zero 
and  maintaining  a  x  pulse  width,  s(t)  becomes: 

s(t)  =  cos{a.0t  +  Df  [.43t  +  ^  sin(^)  +  ^  sin(^)]} 

.  rect  [— t^4-] 

A  plot  of  s(t),  with  a  cosine  window  function  applied  to  it,  is 
shown  in  Figure  13. 

One  advantage  this  new  non-linear  FM  has  over  previously  examined 
non-linear  FM  signals  is  that  its  spectrum  can  be  derived  from  a 
closed-form  analysis.  The  following  is  a  complete  derivation  of  this 


spectrum. 


H3 . 333 
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Figure  13.  Plot  of  New  Non-Linear  FM  Received  Signal. 


First,  the  received  signal,  s(t),  with  the  integral  of  the  h3i(t) 
Eigen  as  its  phase  term  and  a  cosine  window  function  included  can  be 
written  as: 


s(t)  =  cos{wQt  +  .43  Dft  + 


.  5T  Df 
2tt 


sin(^-) 


,  . 07T  Pf 
4tt 


si  n  ( 


by  making  the  following  definitions: 
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0.5T  Df 

6l  =  jj)  “ 


and 

0.07T  Df 
^2  =  4tt 

also: 


2tt 

ml  =  “T 


and 


_  4tt 
m2  "  T 


s(t)  can  now  be  rewritten  as: 


s (t)  =  cos {(jQt  +  . 43Dft  +  sin(wmlt)  +  B2  sinU^t)} 


rect( 


t  +  T/4, 
T/2  ' 


Next,  utilizing  Euler's  theorem,  eiJa)t  =  cos  ut  +  sin  cot,  s(t)  can 
again  be  rewritten  as: 

jooct  j0.43Dft  j 3 1  sin(u)mit)  JB2  sin(wm2^)  ,2Tr(t+T/4l\ 

s ( t )  =  Re{e  e  e  e  >  ‘  cos(^-  yI  J-) 
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Expanding  the  last  two  exponentials  in  a  Fourier  series  yields: 


JBj  s1n(<omlt)  _  - 

'  n:  n  e 

n=-eo 


and 


,Jb2  s1,'(“m2t)  .  ”  F  eJ|taWt 


K=-»  " 

where  the  coefficients,  Fn  and  F^,  are  represented  by: 

T 


!  2  JSi  -Jnlu.  t) 

F„  -  T  /  e  e  dt 


T 

2 


and 


Fl,  = 


1 


T 

7 


f  ;  e 

1  T 


j82  cos(u)m2t)  -jK(o)m2t)  dt 
e 


Next,  by  defining  an  and  a2  as  follows: 


f  2rrt 
°1  “ml4  •  ~T~ 


and 


The  equations  for  Fn  and  F^  can  be  rewritten  as: 


F  =  — 
n  2  it 


tr 

f 

-7T 


COS  a 


-  Zctj) 


da. 


and 


fk'st 


IT 

/ 

-ir 


j{g,  cos  a,  -  Za9) 


dot„ 


F  and  F.,  are  now  in  the  form  of  the  easily  recognizable  Bessel 
n  K 

function  of  the  first  kind,  commonly  written  as: 

F„  *  w 

and 

FK  ■  W 

Values  for  the  Bessel  function  exist  in  tabulated  form  (Watson  1966) 
or  can  be  obtained  numerically  (by  computer)  by  evaluating  the  integral 
(Gradshteyn  and  Ryzhik  1980). 

1  77 

J  (B)  =  ^  /  cos(8e  -  Z  sin  9)  d9  [n  -  a  natural  number] 
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The  received  signal,  s(t),  can  now  be  expressed  in  terms  of  the 
Bessel  functions  as: 


s(t)  =  Re  {e 


JWot 


j  0, 


48Dft 


00 

I 

n=-“ 


W 


jnwmjt 


K=-oo 


jk(b2) 


I***1 


} 


.  cos(Mt^jm) .  rect(^i) 

Now,  again  using  Euler's  theorem,  the  real  part  of  s(t)  can  be  written 
as  a  cosine  such  that: 

oo  oo 

s(t)  =  E  Z  Jn(Bx)  Jk(B2)  COS((w0t  +  0. 43Df  +  ncoml  +  Kw^Jt] 
n~~°°  K=-®  ' 

•  cos(^i^i)  •  rectt^.)  ' 


Since  the  Bessel  functions  are  just  constants,  the  application  of 
some  very  simple  Fourier  transform  rules  will  yield  an  expression  for 
the  spectrum,  s(f). 


s(f)  =  l  I  0„(6i)  JKte2)  {Stf  -  (f0  +  0.43Df  +  nfml  +  Kf^)] 
n--°°  K=-°° 

+jJL 

+  5[f  +  (f0  +  0. 43Df  +  nfml  +  Kf^)]}  *  [6(f  -  ^)e  2  +  6(f  +  gf) 

.j  1  -j  ill 

e  2]  e  *  [J  sine  [—)  e  2  ] 
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where  *  denotes  the  convolution,  in  this  case  in  the  frequency  domain. 

Now,  assuming  f  to  be  much  greater  than  the  swept  bandwidth,  Af, 
so  that  the  negative  frequency  contributions  are  negligible  in  the 
positive  region,  and  considering  only  the  magnitude  for  the  positive 
frequencies : 


I  s(f )  | 


OO  00 


E  E 

n=_oo  K  =-<*> 


W 


0K(e2)  {«[f  -  (f0  +  0.43Df  *  nfnl  +  Kfm2)) 


*  *^sinc(S)) 

Using  Fourier  transform  properties,  the  magnitude  of  s(f)  can  be 
expressed  as: 


Js(0| 


oo  oo 

E  £  Ok(B2)  }  sine 

n=-oo  K=  -oo 


(f-f0-^-Q.43Df-nfml-Kfm2)T 
L  2  J 


where  is  a  constant  used  to  determine  the  swept  spectrum  bandwidth, 
Aoj.  Figure  14  shows  a  plot  of  the  magnitude  of  the  spectrum,  s(f), 
obtained  by  utilizing  the  FFT  (fast  Fourier  transform)  and  plotting 
capabilities  of  SAWCAD. 

For  a  matched  filter  pulse-compression  system,  the  matched  filter 
impulse  response  can  be  written  directly  as: 
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h(t)  -  cos  {.0(I  -  t)  +  0.430f  (j  -  t)  +  Sjjp- 

♦  0,021  •  cos(2li^J24l)  .  rect(±_^4) 

A  plot  or  the  matched  filter  impulse  response  is  shown  in  Figure  15. 

The  desired  output,  the  auto-correlation  function  g(t),  can  be 
obtained  by  several  different  methods.  The  first,  but  most  difficult, 
method  is  to  actually  perform  the  convolution  and  obtain  the  auto¬ 
correlation  function  directly.  In  which  case: 

g(t)  =  s(t)  *  h(t) 

Since  for  a  matched  filter, 


h(t)  =  s ( t  -  t) 


then 


g(t)  “  s ( t)  *  s (t  -  t) 

A  second  method  of  obtaining  the  desired  output  is  to  multiply 
the  signal  spectrum  and  the  matched  filter  spectrum  to  obtain  an 
output  spectrum,  G(w).  Then,  the  inverse  Fourier  transform  will  yield 
the  desired  output,  g(t).  In  equation  form,  it  looks  like: 


G(w)  =  S(w)  ’  H(w) 
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and  then 

g(t)  =  r1  {G(o))} 

This  process  was  accomplished  on  SAWCAD,  and  G(w)  is  shown  plotted  in 
Figure  16  and  g(t)  in  Figure  17.  A  quick  comparison  of  this  auto¬ 
correlation  function  with  that  of  the  linear  FM  will  show  the  reduction 
in  si  delobe  levels.  A  complete  and  thorough  comparison  will  be 
accomplished  in  a  later  section  of  this  paper. 


Comparison  of  Properties  with  Linear  and 
Non-Linear  FM  Systems 

In  any  radar  system,  there  are,  in  general,  three  performance 
parameters  that  are  of  prime  importance.  They  are: 

1.  Bandwidth  -  Range  resolution  is  determined  by  the  frequency 
structure  or  spectrum  of  the  signal.  For  a  given  spectral 
shape,  range  resolution  is  proportional  to  1/F  where  the 
bandwidth  of  the  signal  is  F  cycles. 

2.  Time  Duration  -  Velocity  (that  is,  the  radial  component  of 
target  velocity)  resolution  is  determined  by  the  time  structure 
or  envelope  of  the  signal.  For  a  given  envelope  shape, 
velocity  resolution  is  proportional  to  1/T  cycles,  where  the 
time  duration  of  the  signal  is  T  seconds.  Finest  velocity 
resolution  is  obtained  with  a  signal  which  has  a  rectangular 
envelope. 

3.  Received  Signal,  q(t).  Energy  -  Target  detectability  is 
determined  by  the  ratio  of  received  signal  energy  to 
receiver  noise  power,  S/N  ra*io. 

Of  course,  all  three  of  these  parameters  are  related  to  each  other 
through  time-frequency  domain  transforms. 

In  this  section,  primary  emphasis  will  be  given  to  comparing  the 
properties  of  the  matched  filter  output  (auto-correlation  function) 


FREQUENCY  In  (-Hz 


Figure  17.  New  Non-Linear  FM  Auto-Correlation. 
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of  linear  and  new  non-linear  FM  pulse-compression  systems.  The 
auto-correlation  function  can  be  thought  of  as  a  summary  of  the  three 
parameters  listed  previously.  In  all  data,  it  is  assumed  that  the 
input  to  the  matched  filter,  r(t),  is  ideal,  and  contains  no  Doppler 
frequency  shifts. 

In  this  first  comparison,  the  auto-correlation  function  of  the 
new  non-linear  FM,  Figure  18,  is  compared  to  the  auto-correlation 
function  of  linear  FM,  Figure  19.  Table  2  lists  the  characteristics 
of  each  system  for  a  time-bandwidth  product  of  4.8. 


TABLE  2 

CHARACTERISTICS  OF  NEW  NON-LINEAR  FM 
AND  LINEAR  FM  WITH  A  TIME-BANDWIDTH  PRODUCT  OF  4.8 


NEW  NON-LINEAR  FM 

1  LINEAR  FM 

■  - 

Center  Frequency,  fQ  (MHz) 

70 

70 

Pulse  Length,  t  (micro  sec) 

.42857 

.42857 

Bandwidth,  Af  (MHz) 

11.14 

11.14 

Time-Bandwidth  Product,  TAf 

4.8 

4.8 

Time  Sidelobe  Level  (dB) 

-12.623 

-23.5 

(3  dB)  Compressed  Pulse  Width,  T 
(micro  sec) 

.055 

.078 

Actual  Compression  Rate,  t/T 

7.8 

5.5 

rnocr-i*Nir--o3»  ncHS»rn» 
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A  brief  review  of  the  characteristics  listed  in  Table  2  shows 
that  the  new  non-linear  FM  has  higher  sidelobes,  -12.6  dB,  than  the 
linear  FM,  -23  dB,  It  is  obvious  that  at  this  particular  time- 
bandwidth  product,  the  linear  FM  would  still  be  better.  Also  of 
interest  is  the  fact  that  the  new  non-linear  FM  achieved  an  actual 
compression  ratio  of  7.8  as  compared  to  5.5  for  the  linear  FM.  This 
could  possibly  be  exploited  to  yield  higher  range  resolution. 

In  this  next  comparison,  a  cosine  window  function  has  been  added 
to  both  the  new  non-linear.  Figure  20,  and  the  linear  FM,  Figure  21, 
signals  and  matched  filter.  Table  3  summarizes  the  characteristics 
of  each  system  for  a  time-bandwidth  product  of  4.8. 

TABLE  3 

CHARACTERISTICS  OF  NEW  NON-LINEAR  FM  AND  LINEAR  FM  WITH 
A  TIME-BANDWIDTH  PRODUCT  OF  4.8  AND  A  COSINE  WINDOW  ON  ALL  SIGNALS 


NEW  NON-LINEAR  FM 

LINEAR  FM 

Center  Frequency,  f  (MHz) 

70 

70 

Pulse  Length,  t  (micro  sec) 

,42857 

.42857 

Bandwidth,  Af  (MHz) 

11.14 

11.14 

Time-Bandwidth  Product,  iAf 

4.8 

4.B 

Time  Sidelobe  Level  (dB) 

:  -20 

:  -17.5 

(3  dB)  Compressed  Pulse  Width,  T 
(micro  sec) 

.08385 

.11310 

Actual  Compression  Ratio,  t/T 

5.11 

3.79 

Figure  20.  Plot  of  Auto-Correlation  Function  for  New  Non-Linear  FM 
with  Cosine  Window,  T^f  s  4.8. 


TIME  In  mlcroSECOf'Os 


Figure  21.  Plot  of  Auto-Correlation  Function  for  linear  FM  with 
Cosine  Window,  iAf  =  4.8. 
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The  addition  of  the  cosine  window  has  the  effect  of  reducing  the 
auto-correlation  time  sidelobes,  but  at  the  same  time  increasing  the 
compressed  pulse  width.  The  window  has  caused  these  sidelobes  to  be 
masked  and,  therefore,  not  detectable.  A  comparison  of  the 
characteristics  listed  in  Table  3  shows  that  the  new  non-linear  FM 
now  has  lower  sidelobes,  -20  dB,  than  the  linear  FM,  -17.5  dB.  Also, 
the  new  non-linear  FM  achieved  an  actual  compression  ratio  of  5.11  as 
compared  to  3.8  for  the  linear  FM.  In  this  case,  the  new  non-linear 
FM  would  perform  better  in  a  situation  where  range  resolution  is  of 
primary  concern. 


Design  Procedures  and  Considerations 
The  complete  design  of  any  surface  wave  device  has  to  take  into 
account  the  input  signal,  in  both  the  time  and  frequency  domains,  and 
also  the  various  se..ond-order  effe  s  in  order  to  satisfy  demanding 
system  requirements.  The  general  design  considerations  for  dispersive 
devices  (both  linear  and  non-linear  FM  are  dispersive)  has  been 
examined  (Ristic  1983)  and  will  he  reviewed  here. 

The  first  step  in  designing  a  sampled  SAW  device  is  to  determine 
the  electrode  (finger)  placement.  For  dispersive  devices,  this 
placement  must  be  such  that  the  waveform  is  real.  This  can  be 
accomplished  by  oetermining  when  the  signal  phase  is  either  C  or  it,  as 
illustrated  in  Figure  22  (Ristic  1983). 
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Figure  22.  Determining  Electrode  (Finger)  Placement  at  Time  tm- 


The  sample  time,  tm>  is  determined  by  solving  the  equation:  sin  <f>  (t) 
=  0.  As  an  example,  for  linear  FM: 

♦  (t)  •  2ir  (f0  t  +  t2) 

Mow  solving  sin  <p  (t)  =  0  yields: 

<p  (tm)  =  mrr  +  c 

where  m  is  the  sample  number  and  c  is  an  unknown  constant.  The 
constant  c  is  determined  to  put  the  time  t  at  the  origin  of  the 
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time  scale,  t  =0.  If  M  is  the  total  number  of  samples  (electrodes) 
desired,  then: 

*  (tm}  =  271  (fo  tm  +  If  tm2).=  ^  +  c 

and 

n  if  +  ,  Af  .  2\  _  M 
2  ^fo  ^  2T  ^  m  "  I 

Therefore: 

M  =  2  f0  T  +  1 

Now,  solving  for  t  : 


Af  ,M 


M  +  (i  .  (n.  m)1^} 

f.  T  1 


This  now  determines  the  time  when  the  linear  FM  signal  should  be 
sampled  to  obtain  only  the  real  amplitude  peaks.  This  is  a  form  of 
non-uniform  sampling,  and  is  accomplished  in  the  following  manner. 
The  signal  waveform  is  sampled  according  to: 


h(t)  =  cos  <*>  (t)  i[6  cp( t)  -  mir] 
m 

where  <j>  (t)  a  ^t  +  |y  t2. 


In  order  to  account  for  the  fact  that  the  major  contribution 
to  the  spectrum  occurs  at  the  lower  frequencies,  the  device  will  need 
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to  be  amplitude  apodized  inversely  witw  respect  to  the  instantaneous 
frequency. 

In  equation  form,  this  looks  like: 


6(t  -  t  ) 

h'(tl  =  cos  * (tl  i  trrgr 


where  (tm)  is  the  instantaneous  frequency  at  sample  time  tm>  This 
can  be  rewritten  as: 


h'(t)  =  I  cos  $  (tm) 
m 


<5(t  -  tj 


This  equation  is  now  in  the  form  of  a  general  relationship  that 
can  be  used  to  determine  the  sample  position  and  amplitude  for  any 
dispersive  device.  It  should  be  noted  that  the  amplitude  apodization 
stated  above  would  yield  a  constant  amplitude  signal  under  conditions 
in  which  no  other  losses  are  present.  However,  it  has  been  noted 
(Ristic  1983)  that  losses  which  increase  with  frequency  exist  and  are 
predominant  in  practical  devices.  Therefore,  the  required  amplitude 
apodization  will  be  an  increasing  function  of  frequency. 

In  applying  the  above  relationship  to  the  non-linear  FM 
introduced  in  this  paper,  it  is  noted  that  the  equation  for  tm  is 
quite  complex.  As  an  alternative  to  the  equation,  the  sample  times 
(t  )  can  be  located  by  detecting  the  zero-crossing  times  and  then 
sampling  at  the  midpoint  between  each  pair  of  zero-crossing  times. 
This  can  be  accomplished  by  sampling  the  signal  at  a  very  high  rate 
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to  locate  the  zero-crossings.  The  above  procedure  will  establish 
sampling  times  and  amplitudes.  It  should  also  be  noted  that  when 
the  electrodes  are  laid  down  on  the  substrate*  it  is  desirable  to 
maintain  a  50%  duty  cycle.  That  is,  to  maintain  equal  electrode  and 
spacing  widths.  This  will  require  that  the  finger  width  be  reduced 
proportionally  to  the  increase  in  frequency.  A  sketch  depicting  this 
effect  is  shown  in  Figure  23. 


Figure  23.  Sketch  of  Electrode  (Finger)  Placement  for  New 
Non-Linear  FM. 
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Some  second-order  effects  have  been  considered  in  the  design  of 
a  dispersive  SAW  device  (Jones  et  al.  1972}  and  will  be  reviewed 
here.  Basically,  choice  of  substrate,  beam  spreading,  finger  placement 
accuracy,  weighting  to  reduce  sidelobes,  cascading,  multiple  transit 
reflections,  temperature  effects  and  input  admittance  and  matching 
must  all  be  considered  when  actually  building  the  device.  The 
majority  of  these  secondary  considerations  have  already  been 
incorporated  into  SAWCAD. 

In  an  actual  implementation  of  a  dispersive  device,  it  will  be 
necessary  to  include  compensation  to  account  for  the  fact  that,  input 
admittance  varies  with  frequency.  This  compensation  will  be  in  the 
form  of  amplitude  weighting  proportional  to  1/f.  Figure  24  depicts 
the  way  the  transducer  fingers  will  appear.  Before  this  device  can 
be  implemented  using  SAWCAD,  it  will  be  necessary  to  modify  SAWCAD 
STRUCTURE  program  to  enable  it  to  handle  dispersive  devices. 


Figure  24.  Sketch  Showing  Finger  Placement  and  Amplitude  Weighting. 


CHAPTfft  IV 

CONCLUSION  AND  SUMMARY 


Matched  filtering  and  pulse-compression  concepts  have  been 
reviewed  in  the  context  of  radar  signal  processing.  The  application 
of  these  concepts  to  linear  and  non-linear  FM  have  been  examined 
with  the  primary  emphasis  being  given  to  the  auto-correlation  output 
function.  It  was  shown  that  for  most  applications,  linear  FM 
provides  the  best  overall  performance.  One  limitation  that  exists 
in  linear  FM  is  that  the  auto-correlation  sidelobe  level  is  limited 
to  -23  dB  for  high  time-bandwidth  products. 

Non-linear  FM  pulse-compression  is  one  way  to  exceed  this 
sidelobe  limitation.  Several  successful  versions  of  non-linear  FM 
were  reviewed.  The  major  shortcoming  of  these  implementations  is 
that  the  design  procedure  is  iterative. 

A  new  non-linear  FM  pulse-compression  technique  was  introduced 
and  evaluated  for  implementation  on  a  surface  acoustic  wave  (SAW) 
device  using  SAWCAD.  It  was  shown  that  for  some  applications, 
specifically  low  time-bandwidth  product  cases,  this  non-linear  FM 
could  possibly  perform  better  than  linear  FM. 

In  the  final  section,  design  considerations  are  discussed.  The 
major  design  consideration,  finger  (electrode)  placement,  has  been 
discussed.  The  conclusion  being  that  finger  placement  will  have  to 
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be  accomplished  by  determining  the  zero-crossings  and  sampling  at 
the  midpoint  between  respective  zero-crossings.  It  was  also  shown 
that  amplitude  weighting  will  be  necessary  to  counter  the  effects 
of  increased  admittance  with  increasing  frequencies. 

In  final  conclusion,  the  implementation  of  this  new  non-linear 
FM  using  SAWCAD  will  require  that  SAWCAD  be  modified  to  handle  any 
dispersive  device  design. 
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